Hybrid exchange-correlation functional for accurate prediction of the electronic and 

structural properties of ferroelectric oxides 
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Using a linear combination of atomic orbitals approach, we report a systematic comparison of 
various Density Functional Theory (DFT) and hybrid exchange-correlation functionals for the pre- 
diction of the electronic and structural properties of prototypical ferroelectric oxides. It is found 
that none of the available functionals is able to provide, at the same time, accurate electronic and 
structural properties of the cubic and tetragonal phases of BaTiOa and PbTiOs. Some, although not 
all, usual DFT functionals predict the structure with acceptable accuracy, but always underestimate 
the electronic band gaps. Conversely, common hybrid functionals yield an improved description of 
the band gaps, but overestimate the volume and atomic distortions associated to ferroelectricity, 
giving rise to an unacceptably large c/a ratio for the tetragonal phases of both compounds. This 
super-tetragonality is found to be induced mainly by the exchange energy corresponding to the 
Generalized Gradient Approximation (GGA) and, to a lesser extent, by the exact exchange term of 
the hybrid functional. We thus propose an alternative functional that mixes exact exchange with 
the recently proposed GGA of Wu and Cohen [Phys. Rev. B 73, 235116 (2006)] which, for solids, 
improves over the treatment of exchange of the most usual GGA's. The new functional renders an 
accurate description of both the structural and electronic properties of typical ferroelectric oxides. 

PACS numbers: 71.15.Mb, 71.15.Ap, 71.15.Nc, 77.84.Dy 



I. INTRODUCTION 

Ferroelectric (FE) oxides constitute an important class 
of multifunctional compounds, attractive for various 
technical applications in fields such as electronics, elec- 
tromechanical energy conversion, non-linear optics or 
non- volatile data storage. P, Q Since the early 1990's, 
these compounds were the subject of numerous first- 
principles studies based on Density Functional Theory 
(DFT), most of which were performed within the usual 
Local Density Approximation (LDA) and, to a lesser ex- 
tent, the Generalized Gradient Approximation (GGA). 3, 
0] Although many successes have been achieved, the typ- 
ical inaccuracies of these popular approximations impose 
considerable limitations. 

The LDA is well known to underestimate the lattice 
constants. While the typical underestimation of only 1 
or 2 % of the experimental value is in many contexts re- 
garded as acceptable, it reveals itself quite problematic 
in the study of FE oxides whose structural instabilities 
are extremely sensitive to volume changes and, usually, 
suppressed significantly for the LDA equilibrium lattice 
constants. 3 In some cases, calculations are performed at 
the experimental volume to obviate this problem. How- 
ever, such a method is not fully consistent and is un- 
feasible when experimental data are unavailable. Also, 
in the study of epitaxial multilayers in which different 
materials alternate, because the LDA errors pertaining 
to the lattice constants of the different compounds may 
be slightly different, it is impossible to impose a correct 



epitaxial strain simultaneously to all the layers. 

In most cases, the GGA constitutes a significant im- 
provement over the LDA, although it has a tendency to 
overcorrect the LDA error for the equilibrium volume, 
thus leading to overestimations. In the case of ferro- 
electrics, however, the well known GGA functional of 
Perdew, Burke and Ernzerhof (GGA-PBE) f?\ performs 
significantly worse than the LDA, yielding a wrong super- 
tetragonal structure in BaTiOa and PbTiOa.Q It was 
only very recently that Wu and Cohen [§] proposed a 
modified GGA functional (GGA-WC) that is accurate 
for BaTiOs and PbTiOa and thus opens the door to a 
straightforward and reliable description of the structural 
properties of FE oxides. The so-called Weighted Density 
Approximation (WDA) [T^ has also been shown to con- 
stitute an improvement over the LDA.[1, |Tl[ However, 
since this approach is not available in the most widely 
used DFT codes, it has been applied only marginally. 

It is also well known that the LDA and GGA usually 
lead to a significant underestimation of the band gaps 
of the Kohn-Sham electronic band structure (often by 
as much as a 50 %). 12] This does not a priori consti- 
tute an intrinsic failure: Kohn-Sham particles are ficti- 
tious independent particles resulting from a mathemat- 
ical construction, with no other formal link to the real 
interacting electrons of the system than the requirement 
to reproduce the same total density, and no guarantee to 
exhibit the same energy spectrum. However, although 
"exact" Kohn-Sham DFT should internally correct for 
the so-called band-gap problem in order to provide accu- 
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rate ground-state properties, that is not necessarily the 
case for approximate functionals. For example, the op- 
tical dielectric constant of insulators is usually badly re- 
produced within the LDA and GGA due to the absence 
of ultra non-local dependence of the exchange-correlation 
kernel. [H [13 The band-gap problem becomes even more 
critical in, for example, the emerging field of magnetic 
ferroelectrics, where the LDA and GGA often render er- 
roneous predictions of metallicity for systems that are ac- 
tually FE insulators [H, [H, [13] , thus precluding further 
investigations. Also, pathological situations may occur in 
the study of metal/FE interfaces: Because of the band- 
gap underestimation, the LDA often predicts the Fermi 
level of the metal to be erroneously aligned with the con- 
duction bands of the FE, instead of being located within 
the gap. As a consequence, the LDA renders an unre- 
alistic charge transfer from the metal to the FE, which 
in turn results in incorrect predictions for the properties 
of the interface or FE thin film under study. In view 
of all these failures, there is clearly a need for improved 
alternative functionals. 

The so-called "hybrid" functionals jl^ that com- 
bine Hartree-Fock (HE) and DET, such as B3PW [11], 
B3LYP [13] and Bl, [2l|,[22] are very popular in Quantum 
Chemistry as they have been shown to provide accurate 
description of the atomization energies, bond lengths, 
and vibrational frequencies, together with good energy 
spectra for most molecules. (H IH, [H, [^, [H, ^ In 
the last decade, these hybrid functionals have been also 
increasingly applied to solids. [13, [Ml, [H, [H, [H, [H, [s^, 
[btI . [ssI ] As it was for instance illustrated on Si, B3LYP 
constitutes an interesting alternative to LDA and GGA, 
as it provides excitation energies in much better agree- 
ment with the experiment. [3^ Further, B3LYP has also 
been shown to significantly improve over the LDA re- 
sults for magnetic oxides. ^33, j39, 40, 41, 4^] Such hybrid 
functionals have been recently applied to ferroelectrics 
and, indeed, they provide a much better description of 
the electronic band gaps. [13] However, a more exhaus- 
tive analysis that includes investigation of the structural 
and FE properties is still missing. 

In this paper, we consider two prototypical ferro- 
electrics, barium titanate (BaTiOa) and lead titanate 
(PbTiOa), and show that, although B3LYP and Bl give 
an acceptable description of the electronic properties, 
they fail to reproduce the correct tetragonal polar phase 
of these two compounds due to the super-tetragonality 
introduced by the GGA exchange part and, to a lesser 
extent, the exact exchange contribution. Then, we resort 
to the GGA recently introduced by Wu and Cohen, which 
is known to be accurate for the structural properties of 
ferroelectrics, and construct an alternative hybrid func- 
tional (Bl-WC) that, while remaining accurate for the 
structural properties, is also accurate in the prediction of 
the band gap and electronic structure. This functional 
opens the door to an improved theoretical description of 
both the electronic and structural properties of this im- 
portant class of materials, as required in the especially 



challenging cases mentioned above. 

The outline of the paper is as follows. In Section II 
we briefly introduce the concept of hybrid functional. In 
Section III we describe the technical details of our cal- 
culations. In Section IV we report the results obtained 
using the available functionals (LDA, GGA, B3LYP, Bl 
and HE), for the different phases of BaTiOa and PbTiOg 
and show that none of them gives simultaneously good 
structural and electronic properties. In the second part of 
the paper, we construct a different Bl-WC hybrid func- 
tional by mixing the exact exchange with the GGA-WC 
functional. This functional is shown to provide an ac- 
curate description of both the electronic and structural 
properties of the ferroelectrics considered. 



II. HYBRID FUNCTIONALS 

DFT is a priori an exact theory but relies, within 
the Kohn-Sham formalism, on an universal exchange- 
correlation energy functional, Exci whose explicit form 
is unknown and that must be approximated. In 
the search of accurate exchange-correlation Kohn-Sham 
functionals, an additional difficulty arises from the fact 
that Excj obtained from the integration of the exchange- 
correlation hole, is not equal to the exchange-correlation 
energy of the many-electron interacting system. The dif- 
ference originates in the transfer of part of the many- 
body kinetic energy Txc to the exchange-correlation 
term within the Kohn-Sham formalism. This additional 
kinetic contribution makes Exc even more complex to 
estimate that it could be a priori expected. 

An exact relationship is however satisfied by Exc, 
known as the adiabatic connection formula (or coupling- 
strength integration formula) . [44*] Let us consider the fol- 
lowing family of Hamiltonians with different electron- 
electron interactions governed by a single parameter A 
varying from (nonintcracting limit) to 1 (fully interact- 
ing limit): 

H,iiX) ^ t, + XUee + VX (1) 

where Te is the kinetic energy operator, Uee is the 
electron-electron potential energy operator and vx is de- 
fined in such a way that all these Hamiltonians produce 
the same ground-state density n. For A = 1, we have the 
fully interacting system with vx=i — Vext, the external 
potential, while for A = 0, we have the non- interacting 
system with vx=o = vkSi the Kohn-Sham potential asso- 
ciated to the density n. It can be shown that the Kohn- 
Sham exchange-correlation energy corresponds to the av- 
erage of the exchange-correlation potential energy for A 
ranging from to 1, the integration over A generating 
the kinetic part of Exc, which is known as the adiabatic 
connection formula [3 [IH] : 

Exc^ I dXExcx- (2) 
Jo 
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FIG. 1: Kohn-Sham exchange-correlation energy Exc (light 
gray area ) is obtained as the potential energy of exchange- 
correlation of the truly interacting system Exc,\^i modified 
by a positive quantity of the transfered kinetic energy Txc 
(dark gray area) all along the path of integration of the cou- 
pling constant A. 



Exc.x is the potential energy of exchange-correlation at 
intermediate coupling strength A, 



EXC,X = {^X I Uee \ ^x) - Eh 



(3) 



where Eh is the classical electrostatic Hartree energy and 
Tpx are the ground-state wave functions at the coupling 
strength A. 

This result is illustrated in Fig. [TJ At A = 0, from Eq. 
H Exc,x=o = Ex where Ex = {ipo I Uee \ V'o) - Eh 
is the exact exchange energy of the system defined as 
within the Hartree-Fock method but on the basis of the 
Kohn-Sham wavefunctions ipo- At A = 1, Exc,x=i in 
Eq. [3] corresponds to the many-body potential energy 
of exchange-correlation of the fully interacting system. 
The Kohn-Sham exchange-correlation energy, Exc, is 
provided by Eq. [2] and corresponds to the light gray 
area in Fig. [T] It differs from the many-body exchange- 
correlation energy of the fully interacting system by T^c 
(dark gray area in Fig.[T]) corresponding to the transfer of 
many-body kinetic energy along the path of integration 
over A between and 1. 

From Fig. [T] it appears that Exc could be obtained 
from mixing of exact exchange energy [Ex — Exc,x=o) 
and many-body {Exc,x=i) exchange-correlation poten- 
tial energy (end-points of the integration path over A). 
The exact value of this mixing is however unknown and 
depends on the shape of the curve describing the evo- 
lution of Exc,x with A. If the evolution was perfectly 
linear, we would simply get the exact result, 



Exc — 2^Ex + Exc,x=i) 



(4) 



The simple half-and-half hybrid functional pro posed 
by Becke is based on Eq. 0] and it is defined as [iSl l46j : 
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More sophisticated mixing schemes were further pro- 
posed to better capture the evolution of Exc,x with A. 

B3PW and B3LYP are based on the following mixing 
scheme : 
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where A=0.2, B=0.9, C=0.81 are the three Becke's mix- 
ing parameters which were determined by fitting exper- 
imental data on molecules E^^^ and E^^^ are 
the exchange and correlation energy within LDA, and 
Ef'^ and E^(^^ are the exchange and correlation en- 
ergy within the GGA functional. The B3PW functional 
includes the Becke's GGA exchange j47il ^^"^ ^^i^ GGA 
correlation of Perdew and Wang [4^ , whereas the B3LYP 
functional includes the Becke's GGA exchange \^f\ and 
the GGA correlation of Lee, Yang and Parr 49] for E'^^'^ 
and E'g'^^ respectively. 

The more recent Bl hybrid functional sets B=l and 
C==l simplifying the mixing procedure to the only exact 
exchange mixing parameter A [2l|, : 



Exc - E 



GGA 
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X ) 



(7) 



Fitting the experimental data, Becke determined the val- 
ues of 0.16 and 0.28 for the A parameter depending on 
the choice of the GGA functional used in E^'^^ and 
-£'c'^"^-Ei| Perdew, Ernzerhof and Burke provided a qual- 
itative physical justification for the Bl functional and for 
the empirical value of A parameter. [H^l They showed that 
A«0.25 is to be expected for the atomization energies of 
most molecules, while a larger values may be more ap- 
propriate for total energies of atoms and molecules, and 
smaller values for atomization energies of molecules with 
nearly degenerate ground states of the unperturbed A = 
problem. 



III. TECHNICAL DETAILS 

Our calculations have been performed using the lin- 
ear combination of atomic orbitals (LCAO) method as 
implemented in the CRYSTAL code. 20] Results have 
been obtained using (i) DFT at the LDA and GGA 
levels, (ii) hybrid functionals B3LYP and Bl, and (iii) 
the Hartree-Fock method. For the LDA calculations, we 
used the Dirac-Slater exchange [5l| and the Vosko-Wilk- 
Nusair correlation [52^] functionals. For the GGA calcu- 
lations, we considered the exchange and correlation func- 
tional of Perdew, Burke and Ernzerhof (GGA-PBE) 0] 
as well as that of Wu and Cohen (GGA-WC).[9'] As 
for the hybrids, we considered the conventional B3LYP 
functional and a Bl functional with A = 0.16, using 
Becke's GGA exchange and Perdew- Wang GGA corre- 
lation functionals. pH A different Bl-WC functional that 
we introduce in Section V-B makes use of the Wu-Cohen 
GGA exchange-correlation functional and A = 0.16. 
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We should note here that our calculations for the var- 
ious exchange-correlation functionals differ only in the 
exchange-correlation functional used; that is, all the re- 
maining parameters and approximations, including the 
pseudopotentials, are common to all our calculations. It 
should be noted here that keeping always the same pseu- 
dopotentials, irrespectively of the exchange-correlation 
functional, seems to be the dominant way to proceed in 
the community traditionally working with HF and hybrid 
functionals. The pseudopotentials employed are gener- 
ated from HF calculations, and it can be argued [ijl that 
they constitute a reasonable choice. On the other hand, 
the DFT community working on ferroelectrics has always 
favored an alternative approach that involves the gen- 
eration of the corresponding pseudopotentials for each 
exchange-correlation functional. Being aware of these 
two possibilities, in this work we have performed all our 
calculations following two approaches: (i) we used HF 
pseudopotentials for Ba, Pb and Ti atoms, considering 
explicitly all the electrons of oxygen and (ii) we have 
also performed calculations in which Ti, together with O, 
is treated at the all-electron level. Since, the Ti and O 
atoms play a key role as far as the FE properties are con- 
cerned, such an approach allows us to quantify the mag- 
nitude of the error that might be caused by the use of the 
common HF pseudopotentials. Interestingly, we will see 
that such a pseudopotential approximation is reasonable 
except in the case of the Wu-Cohen GGA functional. The 
HF pseudopotentials and the localized Gaussian-type ba- 
sis sets used here are those from Ref. [43] . The basis sets 
include polarization d-orbitals for O ions and were op- 
timized for both BaTiOa and PbTiOs. For the calcula- 
tions of the cubic phases of other perovskite compounds 
which we have considered, the HF pseudopotentials and 
basis sets are as follow: Sr (from Ref. [43]), Ca (from 
Ref. [11]), Zr (all electron from Ref. Is^ l, K (all electron 
from Ref. ^55^), and Ta (from Ref. 

Other technical details are as follows. Brillouin zone 
integrations were performed using a 6x6x6 mesh of k- 
points, and the self-consistent-field calculations were con- 
sidered to be converged when the energy changes between 
interactions were smaller than 10~^ Hartree. An extra- 
large predefined pruned grid consisting of 75 radial points 
and 974 angular points was used for the numerical inte- 
gration of charge density. For the cubic phases we have 
optimized the lattice constants, whereas for the tetrag- 
onal phases we have performed, except where other- 
wise indicated, full optimizations of the lattice constants 
and atomic positions. The optimization convergence of 
3x10^^ Hartree/Bohr in the root-mean square values of 
forces and 1.2x10^^ Bohr in the root-mean square val- 
ues of atomic displacements was simultaneously achieved 
for both atomic position and lattice constant optimiza- 
tions. The level of accuracy in evaluating the Coulomb 
and exchange series is controlled by five parameters. [20j 
The values used in our calculations are 7, 7, 7, 7, and 
14. For computing the spontaneous polarization a denser 
mesh of 10x10x10 fc-points was used. Also, for phonon 




frequency calculations we have used the finite atomic dis- 
placements of 0.0005 A for numerical evaluations of the 
energy derivatives and the energy convergence was in- 
creased to 10~^^ Hartree. 

The quasiparticle band structure in Section IV-A was 
calculated within the GW approximation [53] as im- 
plemented in the ABINIT code.'Hs'] The LDA eigen- 
values were obtained using separable, extended norm- 
conserving pseudopotentials [5^ and used as zero-order 
input in a perturbative fashion. [6^] The LDA calcula- 
tion was performed for a 4x4x4 fc-point mesh inside the 
Brillouin zone with a cutoff energy of 50 Hartree. The 
dielectric matrix Cq^q/ (q) was computed within the ran- 
dom phase approximation (RPA) for 461 G vectors in- 
cluding 200 bands (20 occupied and 180 unoccupied) in 
the summations. The dynamic dependence of e"^ was 
approximated using a plasmon-pole model [6l| fitted to 
the actual calculated values of at two imaginary fre- 
quencies. The self-energy E was obtained by summing 
over 19 special fc-points in the irreducible Brillouin zone, 
and over 400 bands (20 occupied and 380 unoccupied). 
This procedure was repeated iteratively by correcting the 
unoccupied eigenvalues with a scissor operator, match- 
ing the quasiparticle gap from one iteration to construct 
e^^ and the Green's function for the next iteration, until 
convergence of the quasiparticle gap was reached. The 
eigenstates were assumed not to change. The accuracy 
of this procedure is beyond GOWO in which no update is 
done at all and below self-consistent GW in which both 
the eigenvalues and eigenstates are updated. Finally, let 
us note we checked that, in spite of the technical differ- 
ences between the two calculation schemes (e.g. basis 
sets and pseudopotentials), ABINIT and CRYSTAL give 
essentially the same result for the electronic band struc- 
ture of BaTiOa at the LDA level. It is thus meaningful 
to compare the density-functional results obtained with 
CRYSTAL with the GW ABINIT results. 



IV. COMPARISON OF AVAILABLE 
FUNCTIONALS 

BaTiOs and PbTiOa are two prototypical FE 
perovskites. They both adopt at sufficiently high 
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FIG. 3: Indirect E'g and direct Eg band gaps for 



the 



Ti all-electron calculations within different functionals for: 
(a)BaTi03 and (b)PbTi03. are shown in dashed oblique 
lines and Eg in continuous lines. The experimental E^ values 
are also shown. [gH] The pseudo-SIC and GW values for 
BaTiOs are also included. 



temperature a paraelectric Pviim cubic structure as 
shown in Fig. [21 When the temperature is lowered, they 
both undergo a structural phase transition [Tc — 393 
K and 766 K for BaTiOs and PbTiOs, respectively) 
and exhibit a FE PAmm tetragonal phase at room 
temperature. While PbTiOa remains tetragonal down 
to K, BaTiOs undergoes two additional structural 
phase transitions when the temperature is lowered, to 
orthorhombic and rhombohedral phases, respectively. 
The present study focuses on the cubic and tetragonal 
phases of these two representative compounds. 62] In 
this Section, we compare the ability of the available 
DFT (LDA, GGA-PBE, GGA-WC) and hybrid (B3LYP, 
Bl) exchange-correlation functionals to predict correctly 
their electronic, structural and FE properties. 

A. Electronic band structure 

First, we compare the results obtained with the avail- 
able functionals for the electronic band structure. In 
Fig. [31 wc report the indirect {E^) and direct (Eg) band 
gaps for the cubic phases of BaTiOa and PbTiOa, as 



obtained for Ti all-electron calculations. [63[ The values 
were calculated at the corresponding theoretical lattice 
constants. For BaTiOa (resp. PbTiOs) E'^ occurs along 
[R — F] (resp. [X — F]) direction of the cubic Brillouin 
zone (Fig. [21d). For BaTiOa, the results obtained with 
a recently proposed self-interaction-correction scheme 
(pseudo-SIC) [6J] and including GW corrections are also 
shown. The experimental Eg values for BaTiOa [11] and 
PbTiOa [6^ are finally mentioned for comparison. 

On the one hand, it is shown that both LDA and GGA 
(PBE and WC) strongly underestimate the band gaps 
giving rather similar values of Eg and Eg (Fig. [3]). On 
the other hand, HF strongly overestimates these values. 
Only B3LYP and Bl hybrid functionals give irrroroved 
Eg and Eg values as previously reported in Ref. [43| . 

TABLE I: GW corrections to the LDA eigenvalues for the 
bottom of conduction band (AEcb) and the top of valence 
band (AEvb) at different high symmetry points in the cubic 
Brillouin zone. 





F 


X 


M 


R 


AEcsieV) 


1.01 


1.05 


1.03 


1.08 


AEvBieV) 


-0.54 


-0.59 


-0.55 


-0.54 



It is interesting to discuss not only Eg values but also 
the whole dispersion of the electronic band structure ob- 
tained within the different functionals. In Fig. [H we re- 
port the electronic band structure of cubic BaTiOa calcu- 
lated within LDA, B3LYP and including GW corrections. 
Bl-WC results are also shown although they are only dis- 
cussed in Section V-B. The eV energy level corresponds 
to the top of the valence band (VB). 

The group of valence bands in the energy range be- 
tween —5 and eV has a dominant O 2p orbital char- 
acter, whereas the VB levels localized around —11 eV 
are composed mainly of Ba 5p orbitals (see the labels in 
Fig. [4JI). The first five conduction bands (CB's) corre- 
spond essentially to the Ti 3d states and these split in 
two groups of t2g and Cg symmetry, respectively. The ex- 
perimental results of the photoelectron spectroscopy of 
BaTiOa give a VB spread of - 5.5 ± 0.2 eV, Ba 5p or- 
bitals localized around — 12±0.2 eV, and an Eg value of 
^ 3. 3 ±0.2 eV.[7l| The GW results are in good agreement 
with these data, taking into account the experimental un- 
certainties. 

The GW correction to the LDA band structure does 
not simply correspond to a rigid upward shift of the Ti 3c? 
conduction bands. Table [J shows the individual shifts of 
the bottom of CB {AEcb = E^^ - Ej^j^^) and the top 
of VB {AEvB = - E^B^) at various high symme- 

try points of the cubic Brillouin zone. It highlights that 
the increase of the band gap is due both to an upward 
shift of the bottom of CB and a downward shift of the 
top of VB. The Ba 5p VB states are also differently af- 
fected by the GW correction and, as a result, their energy 
reduces with respect to that of the O 2p states. 



6 




Fig. shows the B3LYP band structure, which is in 
much better agreement with the GW band structure than 
the LDA result. In particular, B3LYP and GW render 
a very similar dispersion of the O 2p top VB states. A 
few noticeable differences between the B3LYP and GW 
band structures exist such as (i) a near degeneracy at 
the X point in the CB of B3LYP which is lifted up in the 
case of GW, and (ii) the position of the Ba 5p VB states 
which are higher in energy. 



B. Structural properties 

Next we investigated the structural properties using 
the different functionals. We performed full optimiza- 
tions of the cubic and tetragonal phases. The results of 
these optimizations are shown in Table HIl 

In the cubic phase, the lattice constants are in good 
agreement with those of Ref. [1^. Theoretical results 
from other authors are also shown for comparison. For 
BaTiOa, both the LDA and GGA-WC slightly underesti- 
mate the experimental lattice constant (4.00 A) whereas 
the GGA-PBE, B3LYP, Bl and HF slightly overestimate 
it. For PbTiOs, the value of the cubic lattice constant at 
the FE phase transition (766 K) is 3.969 A. If this value is 
extrapolated to K it reduces to 3.93 A.[7l| Taking this 
extrapolated value as the reference, then again the LDA 
and GGA-WC slightly underestimate the lattice constant 
whereas the GGA-PBE, B3LYP, Bl and HF shghtly over- 
estimate it. The values of the lattice constants obtained 
from Ti all-electron calculations are very close to those 
obtained using Ti pseudopotentials, the largest differ- 
ence (0.006 A) corresponding to the HF case. In sum- 
mary, these optimizations of the cubic phases show that 
LDA produces the largest deviation 1%) compared to 
the experiment for both BaTiOa and PbTiOs while the 
GGA-WC, Bl and HF schemes give all very good values 
for the lattice constants. 

As for tetragonal phase, all the functionals considered 
here give similar atomic distortions when the atomic posi- 



tions are relaxed at fixed experimental lattice constants. 
(Wu et. al. [1| had already observed that this is the 
case for the LDA and GGA functionals.) The atomic 
distortions for all the functionals are comparable with 
the experiment for BaTiOs and slightly overestimated 
for PbTiOg. 

The situation is different when cell parameters and 
atomic positions are relaxed together. The results are 
reported in Table ITTl On the one hand, the LDA under- 
estimates the volumes and atomic displacements as com- 
pared to experiment but gives good values for the tetrag- 
onality c/a. On the other hand the GGA-PBE, B3LYP, 
Bl and HF schemes overestimate the volumes and atomic 
displacements, leading to very large c/a ratios for both 
BaTiOa and PbTiOs. Such a super-tetragonality had 
already been reported for the GGA-PBE functional in 
the case of PbTiOs.fs'l For BaTi03(resp. PbTiOs), the 
largest overestimate in the tetragonality, ~ 6% (rcsp. 
- 19%), is given by the B3LYP (resp. B3LYP and HF). 
The values of the tetragonality given by all functionals 
within the Ti all-electron calculations are only slightly re- 
duced as compared to those obtained from Ti pseudopo- 
tential calculations, with the exception of the GGA-WC 
results for PbTiOa. In that case, the tetragonality is 
strongly reduced from 1.131 to 1.086, a result that comes 
close to the experimental value of 1.071 extrapolated at 
K.[74| Such a special sensitivity of the GGA-WC func- 
tional to the use of HF pseudopotentials deserves further 
analysis. In summary from Table HIl only the GGA-WC 
functional (and, to a lesser extent, the LDA) succeed 
in correctly predicting the right tetragonal structures of 
BaTiOa and PbTiOg. 



C. Ferroelectric properties 

Finally, we also investigated how the different available 
functionals reproduce typical FE properties such as the 
Born effective charges, the spontaneous polarization and 
the FE double-well energy. 



J 
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TABLE 11: Fully relaxed cubic and tetragonal lattice constants a{A), tetragonality c/a, unit cell volume S7o(A"^), and atomic 
displacements dz within different functionals for BaTiOa and PbTiOs. dz relative to Ba(Pb) are given in fractions of c lattice 
constants. The results of the calculations including the Ti pseudopotential are shown in brackets. 











BaTiOs 














PbTiOg 










LDA 


GGA 


HYBRYDS 


HF 


EXP. 


LDA 


GGA 


HYBRYDS 


HF 


EXP 






PBE 


wc 


B3LYP 


Bl 








PBE 


WC 


B3LYP 


Bl 






Cubic phase 




























a 


3.958 


4.035 


3.990 


4.036 


4.008 




A nn 




3.958 


3.923 


3.958 


3.932 


o.yoo 


qq/ 
0. yo 




(3.958) 


(4.UoDj 




( A C\OCi\ 

(4.UoDj 


(A r\r\o\ 
(4.UU0J 


f4 oil) 




(3.896) 


(o.906j 


/o no 1 \ 
(o.yzi ) 


/o K\ 

(o.9oo ) 


(0.929J 


(3.927) 






3.96° 


4.03'' 




4.04'^ 




4.01" 




o.yo 


3.96" 


3.933^ 


3.96" 




o.y± 






3.95T= 


4.UZO 












3.894'' 


o.y i i 














3.937'' 


4.027^^ 












3.890'= 
















3.946"= 




























Tetragonal phi 


ase 


























a 


3.954 


/I m Q 
4.UI0 




o.yyD 


o.yo ( 


3.996 


3.986" 


3.872 


O.oo4 


0.0 ( U 


o.oiy 


o.olz 


3.798 


3,88^ 




(3.951) 


(4.009) 


(3.978) 


(3.988) 


(3.982) 


(3.977) 






(3.828) 


(3.840) 


(3.817) 


(3.805) 


(3.786) 




c/a 


1.006 


1.035 


1.012 


1.066 


1.036 


1.037 


1.010" 


1.041 


1.221 


1.086 


1.277 


1.224 


1.277 


1.071^ 




(1.011) 


(1.043) 


(1.018) 


(1.080) 


(1.045) 


(1.064) 




(1.054) 


(1.236) 


(1.131) 


(1.287) 


(1.238) 


(1.293) 




no 


62.2 


66.9 


63.9 


68.0 


65.7 


66.2 


64.0" 


60.4 


68.8 


62.9 


71.1 


67.8 


70.0 


62.6-^ 




(62.3) 


(67.2) 


(64.1) 


(68.5) 


(66.0) 


(67.0) 




(60.7) 


(69.3) 


(64.1) 


(71.5) 


(68.2) 


(70.2) 




dZTi 


0.011 


0.018 


0.013 


0.019 


0.019 


0.018 


0.015" 


0.037 


0.062 


0.044 


0.076 


0.063 


0.076 


0.040'' 




(0.013) 


(0.019) 


(0.016) 


(0.019) 


(0.020) 


(0.021) 




(0.037) 


(0.064) 


(0.047) 


(0.079) 


(0.065) 


(0.079) 




dzOy 


-0.014 


-0.039 


-0.022 


-0.057 


-0.040 


-0.040 


-0.023" 


0.090 


0.189 


0.121 


0.223 


0.192 


0.225 


0.112'' 




(-0.018) 


(-0.043) 


(-0.025) 


(-0.065) 


(-0.045) 


(-0.056) 




(0.098) 


(0.195) 


(0.141) 


(0.229) 


(0.199) 


(0.233) 




dzo^ 


-0.009 


-0.022 


-0.013 


-0.031 


-0.022 


-0.018 


-0.014" 


0.106 


0.178 


0.133 


0.198 


0.181 


0.199 


0.112'' 




(-0.012) 


(-0.025) 


(-0.016) 


(-0.036) 


(-0.024) 


(-0.027) 




(0.114) 


(0.184) 


(0.148) 


(0.203) 


(0.186) 


(0.204) 





Optimized cubic lattice constants from: "Ref. [lll/Ref. [1], '=Ref. [|], '^Ref. [73],''Ref. [lll.^Cubic and tetragonal lattice 
constants of PbTiOa extrapolated to K from Ref. U^- "Room temperature data from Ref. [gj. ''Room temperature data 

from Ref. [H 



First, we calculated the Born effective charges Z* of 
the optimized cubic phases. The Born effective charges 
were obtained computing the electronic polarization for 
small individual atomic displacements along the z-axis 
with 0|| and atoms located in the z—0 and z=0.5 
planes, respectively (see Fig. [5^). Z* values obtained 
using different functionals are shown in Table IIIII The 
LDA, GGA-PBE, GGA-WC, B3LYP, and Bl render very 
similar results for the effective charges, whereas Z-Jj and 
Zq^^ obtained using the HF scheme are much smaller than 
those given by the other functionals. At the LDA level 
the Z* values are comparable with those of the previous 
pseudopotential calculations. J3J For the different func- 
tionals which we have considered, ZX,-, and ZX have a 
roughly linear dependence with the band gaps Eg, de- 
creasing with Eg values. In comparison to these re- 
sults, we note that the pseudo-SIC method reported in 
Ref. [ill which provides a gap comparable to the experi- 
ment yields surprisingly small values close to HF. For the 
calculations including Ti pseudopotentials, Z^^^p^ and 
Zq^ do not change significantly. There is only a small 
increase in the values of Z^^ and Zq , suggesting that 



the HF Ti pseudopotential will lead to an overestimation 
of the polarization of the FE phases. 

Then, the spontaneous polarization Pg was also cal- 
culated using the Berry phase approach 'tI, [TSj for the 
fully optimized tetragonal phases corresponding to each 
functional (see Table |TTT|. Because GGA-PBE, B3LYP, 
Bl and HF yield an erroneous super-tetragonal struc- 
ture they overestimate the experimental values of P^, 
which are 27 C/m^ for BaTiOg "76] and 0.5-1.0 C/m^ 
for PbTiOa [t^I . The LDA is found to underestimate the 
experimental value of Pg in the case of BaTiOa and to 
produce a more accurate result for PbTiOa. The GGA- 
WC, on the other hand, renders Pg values that are in 
acceptable agreement with experiment for both BaTiOa 
and PbTiOa. The Ps values obtained with Ti pseudopo- 
tentials are typically larger than those obtained including 
all Ti electrons, which is consistent with the overestima- 
tions of the tetragonality and Born effective charges. The 
largest increase for BaTiOs corresponds to the HF results 
(^^120 mC/m^), whereas for PbTiOa the largest increase 
is given by the GGA-WC (-130 niC/m^). 
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TABLE III: Born effective charges Z* calculated at theoretical cubic lattice constants , macroscopic polarization P of the fully 
relaxed tetragonal phases, and the total energy difference AiSt between the optimized cubic and tetragonal phases of BaTiOa 
and PbTiOa. The results including Ti pseudopotentials are shown in brackets. 











BaTiOs 












PbTiOs 










LDA 


GGA 


HYBRYDS 


HF Exp. 


LDA 


GGA 


HYBRYDS 


HF 


Exp. 






PBE 


WC 


B3LYP 


Bl 






PBE 


WC 


B3LYP 


Bl 






Cubic phase 


























7* 


9 7R 


2.73 


2.76 


2.68 


2.72 


2.64 


3.84 


3.87 


3.85 


3.78 


3.83 


3.67 






(2.74) 


(2 71) 


(2 73) 


(2 67) 


(2 70) 


(2.64) 


(3.80) 


(3.85) 


(3.83) 


(3 77) 


(3 81) 


(3.31) 




^Ti(|e|) 


7.11 


7 17 


7 1 n 


u.yo 


u. y4: 


5.74 


6.96 


7 119 


7n 


u.ol 


fi 7Q 


5.82 






(7.29) 


(7.38) 


(7.31) 


(7.15) 


(7.17) 


(5.98) 


(7.15) 


(7.32) 


(7.17) 


(7.01) 


(7.02) 


(6.06) 






-5.65 


-5.73 


-5.66 


-5.35 


-5.55 


-4.41 


-5.77 


-5.83 


-5.81 


-5.58 


-5.60 


-4.43 






(-5.79) 


(-5.91) 


(-5.82) 


(-5.72) 


(-5.75) 


(-4.63) 


(-5.91) 


(-6.02) 


(-5.95) 


(-5.77) 


(-5.81) 


(-4.64) 






-2.12 


-2.08 


-2.10 


-2.04 


-2.05 


-1.99 


-2.51 


-2.53 


-2.52 


-2.50 


-2.51 


-2.53 






(-2.12) 


(-2.10) 


(-2.11) 


(-2.05) 


(-2.06) 


(-1.99) 


(-2.52) 


(-2.53) 


(-2.53) 


(-2.50) 


(-2.51) 


(-2.47) 




Tetragonal 


phase 


























P(C/m^) 


0.20 


0.39 


0.26 


0.48 


0.39 


0.33 0.27" 


0.78 


1.29 


0.98 


1.40 


1.31 


1.39 


0.5-1'' 




(0.24) 


(0.43) 


(0.31) 


(0.53) 


(0.44) 


(0.45) 


(0.85) 


(1.32) 


(1.11) 


(1.43) 


(1.35) 


(1.43) 




l\Et{meV) 


11.6 


40.4 


14.3 


82.1 


48.2 


20.9 


32.6 


248.4 


99.7 


473.7 


284.6 


454.4 






(4.7) 


(49.3) 


(14.5) 


(99.6) 


(59.0) 


(86.2) 


(67.7) 


(292.4) 


(115.5) 


(525.6) 


(331.4) 


(587.7) 





Experimental data from: "Ref. [t^ and 'Ref. [t^]. 



It is interesting to consider the energy difference AE't 
between the paraelectric cubic and FE tetragonal phases, 
as such a quantity strongly correlates with (and can be 
expected to be roughly proportional to) the FE tran- 
sition temperature. The cubic-tetragonal phase tran- 
sition occurs at 393 K (766 K) for BaTi03(PbTi03), 
which corresponds to an energy of 34 meV (66 meV). 
Hence, for BaTiOs the GGA-PBE, B3LYP, and Bl ren- 
der AEt 's that are larger than our reference value, while 
the LDA, GGA-WC and HF schemes give smaller values. 
For PbTiOa, the Ai?t values obtained using the GGA- 



PBE, B3LYP, Bl, and HF schemes are significantly larger 
than what one would expect. These results seem to be 
consistent with the LDA-based first-principles effective- 
hamiltonian studies in the literature, (tsI , which pre- 
dict transition temperatures that are usually lower than 
the experimental ones (by as much as 100 K in the case 
of BaTiOs and PbTiOa). Our results indicate that the 
situation might be improved by using functionals that 
perform better than the LDA in what regards the de- 
scription of the energetics of the system. 



D. Summary 

The results discussed so far clearly suggest that only 
the hybrid functionals can provide an accurate descrip- 
tion of the electronic structure of FE materials, as the 
LDA and GGA typically underestimate the band gap by 
about a factor of two and HF calculations usually ren- 
der large overestimations. At the same time, and un- 
fortunately, these hybrid functionals yield erroneous pre- 
dictions for the atomic structure of the FE phases, as 
best exemplified by the super-tetragonality problem that 
also pertains to the GGA-PBE scheme. In fact, we find 
that only the GGA-WC, and the LDA to a lesser extent, 
provide an accurate description of the structural proper- 
ties of BaTiOs and PbTiOa. As a result, we find that 
none of the available functionals is able to describe with 



acceptable accuracy both the electronic and structural 
properties of these systems. 



V. Bl-WC HYBRID FUNCTIONAL 
A. Origin of B3LYP and Bl failure 

In order to understand the origin of super-tetragonality 
given by the B3LYP functional, we have performed for 
BaTiOa calculations using different values of the three 
Becke's mixing parameters (see Eq. |6]). The results are 
summarized in Table ITVl 

First we note that the percentage of exact exchange {A 
parameter) strongly affects the results of the band gaps 
and c/a ratios, being responsible for the overestimate of 
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TABLE IV; Fully relaxed lattice constants, tetragonality c/a, 
unit cell volume Slo, and indirect band gap Eg [R — T] of 
tetragonal BaTiOs for different values of the Becke's mix- 
ing parameters. The LDA, Bl and Bl-WC (with the exact 
exchange A=0.16 and A=0.25), and experimental values are 
also included. 





A 


B 


C 


a{A) 


c/a 




ElieV) 


LDA 











3.954 


1.006 


62.2 


2.10 




0.2 








3.948 


1.015 


62.4 


3.77 




0.2 


0.9 





3.976 


1.138 


71.5 


3.85 




0.2 





0.81 


3.934 


1.005 


61.2 


3.73 


B3LYP 


0.2 


0.9 


0.81 


3.996 


1.066 


68.0 


3.80 


Bl 


0.16 


1 


1 


3.987 


1.036 


65.7 


3.45 


Bl 


0.25 


1 


1 


3.977 


1.032 


64.9 


4.26 


Bl-WC 


0.16 


1 


1 


3.962 


1.015 


63.2 


3.44 


Bl-WC 


0.25 


1 


1 


3.954 


1.015 


62.7 


4.24 


Exp. [67] 








3.986 


1.01 


64.0 


3.4 [65] 



the tetragonality obtained within HF. Yet, if we fix A 
at the value yielding accurate band gaps {A « 0.2), it 
becomes clear that the main cause for the B3LYP super- 
tetragonality is Becke's GGA exchange contribution {B 
parameter), which also leads to a relatively large value 
for the unit cell volume. Indeed, for B parameter values 
in the interval between 0.0 and 0.9, and no correlation 
contribution (C=0) the tetragonality (resp. unit cell vol- 
ume) increases by ~ 12% (resp. ~ 15%). It is also worth 
to note that the value of B does not affect the calculated 
band gap significantly. Finally, the GGA correlation (C 
parameter) has a small influence on the band gap, tetrag- 
onality and unit cell volume. For C parameter values be- 
tween 0.0 and 0.81, and no exchange contribution {B=0), 
the tetragonality (resp. unit cell volume) decreases only 
by ~ 1% (resp. - 2%). 

This analysis clearly demonstrates that the failures 
of the B3LYP and Bl functionals, that are, the super- 
tetragonality and the severe overestimation of the unit 
cell volume, are essentially caused by Becke's GGA ex- 
change part, and that the tuning of the GGA correlation 
part cannot compensate for such inaccuracies. 



B. Bl-WC hybrid functional 

The above discussion suggests that the super- 
tetragonality problem associated to B3LYP and Bl might 
be overcome by improving for the GGA exchange part. 
Since, as highlighted above, GGA-WC describes the 
structural properties of BaTiOs and PbTiOa quite ac- 
curately, we have built a so called Bl-WC hybrid func- 
tional by mixing exact exchange with GGA-WC follow- 
ing the recent Bl Becke's mixing scheme. This scheme 
only includes the mixing parameter A (see Eq. ^ , with 
values ranging from 0.16 to 0.28 depending of the choice 



of the GGA.[2ll. [2^ Since A monitor the fraction of ex- 
act exchange, the amplitude of the gap is directly linked 
to its value and constitutes good indicator for fixing its 
value. We have performed calculations with A—QA6 and 
A^O.25 for tetragonal phases of BaTiOs (see Table |IV|. 
From Table IIV[ it is clear that going above A=0.16 
percentage produces unrealistically large band gaps for 
BaTiOa. Moreover, Bl-WC with A=0A6 also gives the 
best structural properties. For comparison the calcula- 
tions within Bl functional are also included in Table [TV] 

In order to demonstrate the universality of this choice 
for the class of ABO3 ferroelectrics, we have performed 
full relaxation calculations of the cubic phases for differ- 
ent compounds. These results are shown in Table |Vl 
In all the compounds considered E^(Eg) occurs along 
i? — r(r — r) direction of the cubic Brillouin zone ex- 
cept for PbTiOs and PbZrOa. For PbZrOa, E^(E^) is 
along R - X{X - X), whereas for PbTiOg, E^(E;^) is 
along A" — r(r — F). The agreement with experiment for 
the band gaps and lattice constant is very good for all 
the considered compounds except PbTiOa and PbZrOa. 
For PbTiOs and PbZrOa, we have also performed calcu- 
lations using Durand pseudopotential for Pb.[8^ Chang- 
ing the Pb pseudopotential does not significantly affect 
the values of the band gaps. In the case of PbTiOa, E^ 
increases by 0.2 eV, whereas for PbZrOs, E^ decreases 
by 0.06 eV. Trying to understand the origin of this dis- 
crepancy, it is worth to notice that the band gaps of 
PbTiOa and PbZrOs were not directly measured but ex- 
trapolated to K. The reported values strongly depend 
of the method and the choice of parameters used to ex- 
trapolate the high temperature data. [9^ 

Since Bl-WC with the exact exchange parameter 
yl=0.16 gives good structural and electronic properties 
for the cubic phases of the considered compounds, fur- 
ther, we have calculated the other structural and elec- 
tronic properties of the cubic and tetragonal phases of 
BaTiOa and PbTiOg (see Table Ell- The Bl-WC elec- 
tronic band structure of BaTiOs obtained from the Ti 
all-electron calculation, is shown in Fig. Hj;. The result is 
similar to that of the GW calculation (Fig. [Hi), except 
the position of the Ba 5p states which are higher in en- 
ergy by ^1.5 eV. We notice that even the degeneracy of 
the CB at the X point, which erroneously appeared in 
the B3LYP calculation, is properly lifted up. This de- 
generacy is lifted up at the level of Bl functional, which 
gives very similar band structure results with those of Bl- 
WC. For the Ti all-electron calculations, we obtain struc- 
tural properties in very good agreement with experiment 
for the cubic and tetragonal phases of both BaTiOs and 
PbTiOs. As compared with the GGA-WC results, Bl- 
WC renders cubic lattice constants that are smaller by 
0.5% for both BaTiOs and PbTiOs- For the tetragonal 
phase the unit cell volume is smaller by ~ 1.1% (resp. 
~ 0.8%), while the tetragonality increases by ^ 0.3% 
(resp. ~ 1%) for BaTiOs(resp. PbTiOs). 

Concerning the EE properties, we obtain very good 
values for both BaTiOs and PbTiOs, the Bl-WC Born 
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TABLE V: Fully optimized lattice constants, indirect Eg and 
direct Eg band gaps of different cubic phases within Bl-WC 
hybrid functional. The experimental values are also shown for 
comparison. The results with Ti pseudopotentials are shown 
in brackets for BaTiOs and PbTiOg. 





a 


.(A) 


ElieV) 






Bl-WC 


Exp. 


Bl-WC 


Exp. 


Bl-WC 


Exp. 


BaTiOs 


3.971 


4.00 


3.39 


3.20[65] 


3.45 






(3.97) 




(3 16) 




f3 21) 




PbTiOs 


3.901 


3.93[72] 


2.68 


3.39f93] 


4.23 






(3.90) 




(2.43) 


3.40[66] 


(3.98) 




SrTiOa 


3.880 


3.890" [81J 


3.57 


3.25[82l 


3.91 


3.75f82j 


CaTiOa 


3.834 


3.836'' [83] 


3.59 


3.57[84l 


4.04 




BaZrOa 


4.195 


4.192[85] 


5.24 


5.33[86J 


5.41 




PbZrOa 


4.148 


4.130'= [87] 


3.65 




3.40 


3.83[93J 


SrZrOa 


4.138 


4.109 [88] 


5.51 


5.60[89] 


5.77 




CaZrOa 


4.111 




5.35 




5.65 




KTaOa 


3.971 


3.988'' [90| 


3.41 


3.57[91] 


4.07 


4.40f91J 










3.64[92] 




4.35[92] 



"Extrapolated to K. ''At 600 K. ^Extrapolated to K. 
''Room temperature value. "^Extrapolated to K. 



effective charges being comparable to those obtained with 
the other functionals (except HF) . The energy differences 
AEt between the cubic and tetragonal phases are signifi- 
cantly larger than those obtained from LDA calculations. 
According to the discussion above, in the case of BaTiOs 
such an increase probably implies that a Bl-WC-based 
thermodynamic calculation will render transition tem- 
peratures in significantly better agreement with exper- 
iment than those obtained from the LDA. In the case 
of PbTiOa however, the increase in AEt seems rather 
large, and will probably overcorrect the transition tem- 
peratures. 

As a final check, we also investigated to which extent 
the new Bl-WC hybrid functional performs in predicting 
the vibrational spectra. We have calculated the phonon 
frequencies at the zone center (F point) for cubic phases 
and compared these values with the corresponding val- 
ues within LDA and GGA-WC. Phonon frequency calcu- 
lations in the literature are performed at the optimized 
or experimental lattice constants; thus, for a better com- 
parison, here we report both types of calculations. 

In the cubic perovskite structure at the F point there 
are 12 optical modes: three triply-degenerated modes of 
Flu symmetry and one triply-degenerated silent mode of 
symmetry. Transverse frequencies of these different 
modes are reported in Table IVIIl together with experi- 
mental values in the case of BaTiOa. 94] 

At the optimized lattice constant, for BaTiOa, the best 
agreement with the experimental data is achieved with 
the Bl-WC hybrid functional. It is for the Fi„(r01) fer- 
roelectric soft mode that the dispersion of the results is 
the largest, which is consistent with the well known fact 



TABLE VI; Born effective charges Z* of cubic phases, fully 
optimized tetragonal lattice constants, indirect Eg and di- 
rect Eg band gaps of tetragonal phases, atomic displacements 
dz in fractions of c lattice constant, macroscopic polarization 
P, and energy difference AEt between cubic and tetragonal 
phases for BaTiOs and PbTiOs within Bl-WC hybrid func- 
tional. The results with Ti pseudopotentials are shown in 
brackets. The experimental values are also given. 



BaTiOa PbTiOa 
Bl-WC Exp. Bl-WC Exp. 



Cubic phase 










^Ba{Pb)U^|j 






O.OO 






(2.72) 








^Ti(|e|) 


7.08 




6.81 






(7.11) 








^0||(le|) 


-5.57 




-5.62 






(-5.68) 




(-5.76) 






-2.12 




-2.51 






(-2.08) 




-2.51 




Tetragonal ph 


lase 








1- X \ 
a(A) 


3.962 


3.986 [67| 


3.846 


3.88 [72] 




(3.957) 




(3.810) 




c/a 


1.015 


1.010 \&J} 


1.097 


1.071 [72] 




(1.022) 




(1.154) 






63.2 


64.0 [67] 


62.4 


62.6 [72] 




(63.3) 




(63.9) 




ElieV) 


3.44 


3.40 [65] 


2.83 


3.60[93] 




(3.22) 




(2.66) 






3.73 




4.94 






(3.57) 




(4.76) 




dZTi 


0.015 


0.015 [67] 


0.046 


0.040 [68] 




(0.017) 




(0.050) 






-0.024 


-0.023 [67] 


0.129 


0.112 [68] 




(-0.029) 




(0.154) 






-0.014 


-0.014 [67] 


0.140 


0.112 [68] 




(-0.017) 




(0.158) 




P{C/m') 


0.28 


0.27 [76] 


1.03 


0.5-1.0 [77] 




(0.33) 




(1.19) 




AEtimeV) 


24.0 




110.6 






(24.7) 




(162.0) 





that the FE soft-mode frequency is strongly dependent 
on the unit cell volume. Consequently, the dispersion in 
the computed soft-mode frequencies reduces significantly 
when the calculations are made at the experimental lat- 
tice constant. For PbTiOs, although we cannot compare 
with experimental data, we observe that Bl-WC provides 
results comparable to the other functionals, all of them 
being in rather close agreement. We finally notice that 
the soft-mode (and to a lower extent the Fl„(T02) mode 
of PbTiOs) is strongly affected by the use or not of the 
Ti HF pseudopotential. This further highlights the very 
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TABLE VII: Phonon frequencies (cm~^) at the F point for 
cubic phases of BaTiOs and PbTiOa calculated at optimized 
and experimental lattice constants. The results with Ti pseu- 
dopotentials are shown in brackets. The experimental values 
for BaTiOa are also included. 

BaTiOs PbTiOa 
LDA WC Bl-WC Exp. LDA WC Bl-WC 



a — cLopt 



FiuiTOl) 


75i 


128i 


145i 


soft 


127i 


132i 


146i 




(157i) 


(192i) 


(213i) 




(148i) 


(163i) 


(196i) 


Fiu{T02) 


193 


186 


195 


182° 


145 


141 


138 




(193) 


(186) 


(195) 




(128) 


(115) 


(120) 


Fiu{T03) 


480 


469 


482 


482° 


515 


510 


513 




(474) 


(463) 


(476) 




(509) 


(495) 


(506) 




286 


282 


299 


306'' 


219 


211 


231 




(283) 


(280) 


(298) 




(214) 


(208) 


(229) 



Fiu{T01) 


187i 


155i 


208i 


soft 


139i 


135i 


155i 




(218i) 


(204i) 


(255i) 




(170i) 


(168i) 


(192i) 


Fiu{T02) 


181 


183 


186 


182° 


132 


137 


127 




(180) 


(183) 


(185) 




(107) 


(109) 


(96) 


Fiu{T03) 


456 


462 


468 


482° 


491 


503 


497 




(453) 


(458) 


(463) 




(479) 


(487) 


(481) 


F2u 


281 


281 


295 


306'' 


216 


210 


229 




(278) 


(279) 


(294) 




(213) 


(208) 


(228) 



Experimental data from: °Ref. [ij ''This value was 
measured in tetragonal phase. 



delicate nature of the ferroelectric instability and the ne- 
cessity to be very careful in order to reproduce it prop- 
erly. 

VI. CONCLUSIONS 

In summary, our electronic structure calculations of 
prototypical FE crystals BaTiOa and PbTiOs using the 
most popular exchange-correlation functionals, show that 
it is difficult to obtain simultaneously good accuracy for 



structural and electronic properties. On the one hand, all 
the usual DFT functionals (LDA, GGA-PBE, GGA-WC) 
reproduce the structural properties with various degrees 
of success, the recently introduced GGA-WC being by far 
the most accurate, but severely underestimate the band 
gaps. On the other hand, the B3LYP and Bl hybrid 
functionals correct the band-gap problem, but overesti- 
mate the volumes and atomic distortions giving a super- 
tetragonality comparable to that rendered by the GGA- 
PBE for the tetragonal phases. 

We found that the super-tetragonality inherent to 
B3LYP and Bl calculations is mostly associated to the 
GGA exchange part and, to a lesser extent, to the exact- 
exchange contribution. To bypass this problem, we have 
proposed a different Bl hybrid functional obtained by 
combining the GGA-WC functional with a small percent- 
age of exact exchange (A=0.16). With this Bl-WC, we 
have obtained very good structural, electronic and ferro- 
electric properties as compared to experimental data for 
BaTiOa and PbTiOa. Very good agreement with exper- 
iment is obtained for the lattice constants and the band 
gaps for cubic phases of other perovskite compounds, ex- 
cept the band gaps of PbTiOg and PbZrOg. This differ- 
ent Bl-WC functional thus opens the door to a better 
description of the optical properties of ferroelectrics and 
of metal/FE interfaces for which the tetragonal distor- 
tion of the crystalline cell, the atomic displacements, the 
electronic structure, and in particular the band gap, have 
to be accurately described simultaneously. 
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